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^""i ' Abstract 

An Effective String Rope Model (ESRM) for heavy ion collisions at 
RHIC energies is presented. Our model takes into account baryon recoil 
£N| ' for both target and projectile, arising from the acceleration of partons 

in an effective field, produced in the collision. The typical field strength 
, (string tension) for RHIC energies is about 5 — 12 GeV/fm, what allows 

us to talk about "string ropes". Now we describe an "expanding final 
q h i streaks" scenario @, in contrast to a "non-expanding final streaks" 

^ ■ discussed in Ref. ||. The results show that a QGP forms a tilted 

(]} \ disk, such that the direction of the largest pressure gradient stays in 

1-^ ■ the reaction plane, but deviates from both the beam and the usual 

!> ■ transverse flow directions. The produced initial state can be used as an 

^ , initial condition for further hydrodynamical calculations. Such initial 

$—1 ' conditions lead to the creation of a third flow component. 

The realistic and detailed description of an energetic heavy ion reaction 
requires a Multi Module Model, where the different stages of the reaction 
are each described with a suitable theoretical approach. It is important that 
these Modules are coupled to each other correctly: on the interface, which is a 
three dimensional hypersurface in space-time with normal da 11 , all conservation 
laws should be satisfied, and entropy should not decrease. These matching 
conditions were worked out and studied for the matching at FO hypersurfaces 
in details in Refs. P], Qj. 

1 To appear in the Proceedings of the 30th International Workshop on Gross Properties 
of Nuclei and Nuclear Excitation: Hirschegg 2002: " Ultrarelativistic Heavy Ion Collisions", 
Hirschegg, Austria, 13-19 Jan 2002. 



In energetic collisions of large heavy ions, one-fluid dynamics is a valid and 
good description for the intermediate stages of the reaction. Here, interac- 
tions are strong and frequent, so that other models, (e.g. transport models, 
string models, etc., that assume binary collisions, with free propagation of 
constituents between collisions) have limited validity. On the other hand, the 
initial and final, Freeze Out (FO), stages of the reaction are outside the domain 
of applicability of the fluid dynamical model. For the highest energies achieved 
nowadays at RHIC, hydrodynamic calculations give a good description of the 
observed radial and elliptic flows |£|, £|, ||, ||. 

The initial stages for RHIC energies are the most problematic. Non of 
the theoretical models currently on the physics market can unambiguously 
describe the initial stages (see []7| for a discussion). 



1 Expanding final streaks in ESRM 

In Refs. ||, [lO], § we assumed that the final result of the collision of two streaks, 
after stopping the string's expansion and after its decay, is one streak of the 
length Alf with homogeneous energy density distribution, e/, moving like one 
object with rapidity yj (e/ and yf can be easily found from conservation laws). 
The typical values of the string tension, a, are of the order of 10 GeV/ fm, 
and these may be treated as several parallel strings. We assumed that such 
a final state is due to string-string interactions and string decays, which we 
are not going to describe in our simple model. One of the simplest ways to 
quantitatively take into account string decays is presented in Ref. . 



Now we would like to describe the scenario with expanding final streaks, 
which seems to be more realistic [^]. It is based on the solution for the one- 
dimensional expansion of the finite streak into the vacuum, which can be found 
[0] generalizing the description in 

The initial condition is 
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where _R_, R + are the borders of the system. 

For the forward-going rarefaction wave, v > 0, generated on the right end 
of the initial streak, we infer that the head of this wave (the point where 
the rarefaction of matter starts, i.e., where the energy density e starts to fall 



below eo) travels with the velocity (tanhy — c o)/(l — cotanhy ) to the left. 
On the other hand, the base of the rarefaction wave (the point where e starts 
to acquire non- vanishing values) travels with light velocity v — 1 to the right. 
The solution for a simple wave has a similarity form []12| |, i.e., the profile of the 
rarefaction wave does not change with time when plotted as a function of the 
similarity variable (± = 2= ^~ f ^ ± . Thus, 
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Similarly for the backward-going rarefaction wave, v < 0, we get 
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This simple analytical solution is valid as long as the two rarefaction waves 
did not overlap in the middle of the system. Further evolution becomes more 
complicated and doesn't have a similarity form any more. 

Thus, based on the solution presented above we can have a more advanced 
description of the final streaks. Let us assume that the homogeneous final 
streaks, with some e/, yf, are formed in the Center-of-Rapidity frame (CRF) 
when the larger of initial streaks reach the rapidity yf. Up to this point the 
fluid cell trajectories are the same as those for non-expanding scenario ||, but 
after the homogeneous final streak is formed, it starts to expand according to 
the simple rarefaction wave solution, eqs. (|3|-§]) [0- Figure [T] (B) shows new 
trajectories of the streak ends (compare to (A) from Ref. ||), and Figures [1] 
(C,D) present energy density and rapidity profile for this expanding streak. 

Such an initial state with expanding streaks will also help to avoid the 
problem which may cause the development of numerical artifacts, namely a 
step-like in the beam direction initial energy density distribution (output of 
ESRM): it has a jump from e/ = const inside the matter, to in the outside 



vacuum (of course, where there is also a jump of e as a function of x, y, but 
it is much smoother and this is not the direction of the initial expansion). In 
order to avoid (or at least to suppress) this effect, it was proposed in Ref. |H| 
to smooth over initial energy density distribution, for example by a Gaussian 
shape. Our simple analytic solution smoothes out this jump in a natural way. 

Fig. |2] shows the energy density distribution in the reaction plane for the 
RHIC collisions calculated in the scenario with expanding final streaks. The 
QGP forms a tilted disk for b ^ 0. Thus, the direction of fastest expansion, 
the same as the largest pressure gradient, will be in the reaction plane, but 
will deviate from both the beam axis and the usual transverse flow direction. 
So, the new "third flow component" |l4j] appears in addition to the usual 



transverse flow component in the reaction plane. For non-expanding streaks 



this was shown in ||15|| . If we let final streaks expand, this smoothes out the 
picture, but most of the matter, nevertheless, keeps a similar energy density 
profile and velocity distribution. 
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Figure 1: A) The typical trajectory of the ends of two initial streaks corre- 
sponding to the numbers of nucleons n\ and n 2 , e = 65 GeV/nucl, A = 0.09 
(see ESRM description, for example ||). Stars denote the points where 
2/1,2 = Vf- From t = to until these stars, the streak ends move according 
trajectories derived in Ref. || (see eqs. (37, 38)). Then the final streak starts 
to move like a single object with rapidity y/ in CRF. B) The same situation 
as in subplot (A), but for expanding final streak. C) Shows e(z) profiles of 
expanding final streak (th = 5 fm). We can clearly see three regions - two of 
forward and backward rarefaction waves and a middle where the initial energy 
density, ef, is still preserved. D) Shows the rapidity profile of the expanding 
final streak (th = 5 fm). We can clearly see three regions - two of forward and 
backward rarefaction waves and a middle where the initial rapidity, yf, is still 
preserved. 
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Figure 2: Au+Au collision at e = 100 GeV/nucl, {b = 0.5-2 R Au ) , E = T 00 
is presented in the reaction plane as a function of x and z for th = 5 fm/c. 
Subplot A) A = 0.065, subplot B) A = 0.08. The QGP volume has a shape of 
a tilted disk and may produce a third flow component. 



